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Motion About the Stable Libration Points in the 
Linearized, Restricted Three- Body Problem 

Don Mittleman * 


1.0 SUMMARY 

Starting from the well-known differential equa- 
tions of motion for a point particle in the neighbor- 
hood of a triangle libration point (L 4 or L s ) in the 
linearized, restricted problem of three bodies in the 
plane, and after obtaining the finite equations of mo- 
tion, three integral invariants of the motion are 
derived. The Jacobi integral is expressed linearly in 
terms of two of these. The trajectories for varied ini- 
tial conditions are drawn; the motion of a point parti- 
cle, as it moves along its trajectory, can be visualized 
in elementary geometric terms. Each trajectory has 
an exterior bounding curve; its equation is found. An 
approximation to this boundary curve, a well-defined 
ellipse, had been known; the two curves are com- 
pared graphically. For certain initial conditions, 
there is an interior region from which the trajectory 
is excluded; the equation of the boundary of this 
region is found. 


2.0 INTRODUCTION 

The motion of a particle in the neighborhood of a 
triangular libration point in the linearized, restricted 
problem of three bodies in the plane is investigated. 
In section 3, the linearized equations of motion are 
given; in section 4, three invariants of the motion are 
determined; and, in section 5, the trajectories for 
several initial conditions are discussed. All trajecto- 
ries are bounded; thus, in section 6, the equations for 
the exterior bounding curve are derived and, where 
appropriate, the equations for the interior boundary 
curves are given. 


3.0 THE EQUATIONS OF MOTION 

In section 3.1, the differential equations of motion 
for a particle are derived. The Jacobi integral and the 
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Jacobi ellipse are obtained. Generic solutions of the 
differential equations are given in section 3.2 and 
their dependency on initial conditions is shown in 
section 3.3. Various combinations of the constants of 
the system are interrelated; these relationships are 
given in section 3.4. 


3.1 The Differential Equations 

The following discussion parallels closely that 
given in reference 1. Two bodies of finite mass, 
shown as the Earth (£) and the Moon (A/) in the 
following sketch, are assumed to rotate about an axis 
through their center of mass in circular orbits 
lying in a plane tt that is orthogonal to £j. A rotating 
rectangular coordinate system (f^) lies in tt\ the 
Earth and Moon have the fixed coordinates (x j ,0,0) 
and (x 2 , 0,0), respectively, and their center of mass is 
at (0,0,0). The distance between £and A/is chosen to 
be equal to 1; the mass of M is taken as /x 0.5, that 
of £ as 1 — fj.. The choice for the unit of time is such 
that the gravitational constant k 2 = I. 1 



■For the Earth-Moon system, time t * is calculated from this 
dimensionless time t using The equation t* = (T ( l2n)t, where 
T ( is the time for one revolution of the system about the center of 
mass. For T ( = 27.322 days, t* — 4.348 1. 


1 



In the three-dimensional (^t^j) coordinate 
system, the differential equations of motion for a 


point particle are 


d\ 


dVi 

2 ~dT = ^~ 

(1 - *«)($! - *i) M(Si - * 2 ) 

dt 2 

dr 2 

3 r 3 

r l r 2 

d \ 

d 2 W l 

dt i 

2 dT = ^ 1 - 

(1 - M)t7! dVi 

dt 2 

dt 2 + 

3 3 

r l r 2 


rf2f l _ 

(1 - M)fj 


d 2 h 

dt 2 

df 2 

V 

r 2 3 


equations of motion are 


dH 2 drj 2 3 3^3 

— “ 2 dT 4*2 + ^0 - 2 ^2 (la) 


!h + 2 ^2 = 3V3 


df = 4 (1 “ 2M)? 2 + i«2 (IW 


= -r, do 


where z"! 2 = (£ t — x 1 ) 2 + t) 1 2 + ^ 2 andr 2 2 = (^ — 
x 2 ) 2 + t?i 2 + £j 2 . If the potential function 


"(Mi-fi) 'j(«i 2 + ii 2 ) 


t t iL 


Equation (lc), being independent of equations (la) 
and (lb), shows that the motion parallel to the £ 2 
axis is periodic with period 2v. No more will be said 
about this; the remaining discussion concentrates on 
the solutions of the first two equations. 

If equation (la) is multiplied by d| 2 /dr and equa- 
tion (lb) by dqjldt and the results added, the follow- 
ing is obtained. 


is introduced, the equations of motion become 


d2 t i d n i _ a u 

dt 2 dt ~ 9 ?1 

d 2 v i [ 2 £ii = du 

dt 2 dt 9T?! 

d 2 h = u 

dt 2 ^ 1 


If a (£ 2 ,t) 2 ,£ 2 ) coordinate system is defined with 
the origin at L, a triangular libration point, and with 
the axes parallel to the (^t)^^) axes, the linearized 


d$ 2 d 2 | 2 dT) 2 d 2 t? 2 

d^~ Ht 2 


3 d _h Vf 

4 S 2 dt 4 U 


/ d% 2 dn 2 \ 

2 4 2 ^ + ) 


9 dri 2 
+ 4^2-jT 


This equation can be integrated and yields the Jacobi 
integral 



4 2 



2m)I 2 t? 2 



( 2 ) 


2 



where 


where C is a constant to be determined from the ini- 
tial conditions. In particular, three different values 
for Care distinguished: —1, 0, and +1. Since equa- 
tions (la) and (lb) are linear in £ 2 and t) 2 and equa- 
tion (2) is a homogeneous function of £ 2 , tj 2 , d? 2 /t/r, 
and di\i!dt, any other value of C can be reduced to 
one of these three by properly scaling the (£ 2 ,t) 2 ) 
variables. 

3.1.1 Rest trajectories. — Those trajectories for 
which the initial velocity is zero will be called “rest 
trajectories.” For these, C = — 1 and the initial posi- 
tions lie on the conic 

|s 2 2 - 2 m)? 2 tj 2 +f tj 2 2 =\ 

If 0 < fi < 1, the discriminant of the quadratic is 
D »[M(, _ 2m)] 2 - 4(f) < 0 

and the conic is an ellipse which shall be called the 
Jacobi ellipse. 

3.1.2 Simplified form of the equations of motion . — 
To simplify subsequent discussions, an (x,y) coordi- 
nate system is introduced. The origin is taken at L , 
which is also the center of the Jacobi ellipse; the x 
axis coincides in direction with the major axis of the 
ellipse and the y axis with the minor axis. This is 
effectively a rotation of axes with respect to the 
(£ 2 ,t} 2 ) axes through an angle 0, where 9 = 0.5 arctan 
[— V3 (1 — 2/i,)] radians. 

In the (x,j) coordinate system, the linearized 
equations of motion (equations (la) and (lb)) 
become 


x — 2y = Wj 2 x 

(3a) 

y + 2x = co 2 2 v 

(3b) 



and the dots indicate differentiation with respect to 
time. 

In this coordinate system, the Jacobi integral is 
x 2 + y 2 = cjj 2 x 2 + co 2 2 j ' 2 + C (3c) 
and the Jacobi ellipse is 

cj 1 2 x 2 + c o 2 2 y 2 = 1 (3d) 

3.2 Solutions of the Differential Equations 

The solutions to equations (3a) and (3b) are of the 
form 

x = cos \t + K 2 sin Xt 
y = Lj cos Xl + L 2 sin Xt 

Nonzero solutions can be obtained if X satisfies the 
characteristic equation: 



3 



Only the stable solutions are considered; i.e., 
those for which the X are real and not equal to zero. 
Acceptable values of /a are those for which 1 — 

27/la (1 — /a) > 0 or fi < 0.0385207 For the 

Earth-Moon system, /a = 1/82.45 = 0.0121285 . . . 
and the linear theory may be used to describe the 
motion of a satellite near the libration point. 

The solutions to equations (3a) and (3b) are 
among the functions 


3.3 Initial Conditions 

The dependency of K x , K 2 , K 3 , and K 4 on the ini- 
tial conditions follows. If, at time / = 0, x = x g , y = 
y 0 , x = x 0 , and y = y 0 , then 



(5a) 


x = Ky cos Xj t + K 2 sin X^t + A' 3 cos X 2 f + K 4 sin X 2 t 
y = L j cos Xjf + Z , 2 sin X^t + L 3 cos X 2 t + L 4 sin X 2 ? 







CO, 


- X, 


2 


(5b) 


where Ky , . . . , L 4 are arbitrary constants. If these 
trial solutions are to satisfy equations (3a) and (3b), 
the coefficients Ky , . . . , L 4 are related. The solutions 
to equations (3a) and (3b) can be represented by 


K* 


K 2 + h 2 ) x o + 2 y 0 

x 2 2 - *i 2 


(5c) 


x = ^A' 1 cosX 1 t + AjSinXj/j 
+ cos X 2 r + K 4 sin X 2 t ^ 


(4a) 


Ka 



2 h \ 


CO, 



(5d) 


The determination of K t (i = 1, 2, 3, 4) is obtained 
from equations (4a) and (4b) in a straightforward, if 
tedious, calculation. 


y = 



K 2 cos Xyt — Ky sin 



3.4 Some Useful Relationships 


+ M 2 (k 4 cos \ 2 t - K 3 sin X 2 r) (4b) By direct computation, it can be shown that 


where 


and 


My = 


Uy 2 + Xj 2 

2X, 


m 2 = 


“l 2 + ^2 2 

2X n 


coj 2 + co 2 2 = 3 


Xi 2 + x 2 2 = 1 


X1X2 


Mi 


COy 2 + Xy 2 

2X, 


2Xy 

C0 2 2 + Xy 2 


co. 2 + X 2 


2^2 

co- 2 + X- 2 
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and 


Since X,- (/ = 1, 2) is a root of the characteristic 
equation (6), which can be written in the form 


CO 


M.M. =~ 

11 CO- 


Wj 2 + X 2 2X 
2X co 2 2 + X 2 


Alternately, if one starts with the equation 

x - 2y = cjj 2 x (3a) 


the alternate forms for A/j and Af 2 become apparent. 
Finally, from 


2X, 


y + 2x = co 2 j 


(3b) 


"i = 


+ V 


and seeks a solution of the form 

x = /T 1 cos X/ + K 2 sin Xf 
y = L 1 cos Xr + I 2 sin Xt 

the characteristic equation for X is 


X 4 + (o^ 2 + co 2 2 - 4^X 2 + cjj 2 ^ 2 = 0 (6) 


w l w 2 


co 2 2 + Xj 2 


2co^co 2 


X-,c o 


2 2 


1 + 


\\ 2 


2c o. 


co, 


X 2 co 2 


1 + 


From this, it is inferred that 


X 


2 

1 


+ X, 


(coi 2 + co 2 2 - 4) 


or 


Xj 2 + X 2 2 + cjj 2 + co 2 2 = 4 


and, since coj 2 + co 2 2 = 3, Xj 2 + X 2 2 = 1. From 
Xj 2 X 2 2 = Wj 2 w 2 2 , it is inferred that XjX 2 = o)j(u 2 , dis- 
regarding the negative value for the square root. 


2w 1 X 2 




w 2 ( 

> 2 2 + w r 

" lj 

< 2X 2 

co 2 l 

y w i 2 + x 2 

co_i 

l 

w 2 

m 2 

"1 


u 2 
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4.0 INVARIANTS OF THE MOTION 

In section 4.1, two invariants of the motion, Jy and 
J 2 , are derived. A third invariant relationship is given 
in section 4.2, where it is also shown that these three 
invariants can be used to distinguish between 
different trajectories or to identify a “second” trajec- 
tory as a continuation of an earlier identified one. A 
linear relation exists between Jy 2 ,J 2 2 anc ^ Jacobi 

constant; this is derived in section 4.3. 


4.1 The /, and J 2 Invariants 

Suppose a particle P is observed at t = at posi- 
tion (x 0 ,y 0 ) with velocity (x 0 ,y 0 ) . Its trajectory can be 
determined. At some later time fj, the particle Pwill 
be at the point [x(t 1 ,).y(t 1 )] = (xjj^) with velocity 
[-xr(fl) jK*i)] = (-XjJ'i)- However, between t = t^and 
t = ty, you lose track of the particle. Is it possible to 
decide if the particle P observed at time t = fyand the 
particle observed at t = ty are the same? The follow- 
ing discussion partly answers this question. 

Having observed the particle initially at (x 0 ,y 0 ) 
with “initial” velocity (x 0 ,y 0 ), a unique trajectory 
and, in particular, the numbers K 2 , K 2 , and K 4 
are determined. At the later time /j, a particle is ob- 
served at the point with velocity (XjO>j). For 

these “initial” conditions, a unique trajectory and, in 
particular, four numbers AT*, K 2 , Kf, and K 4 are 
determined. The question may then be framed: Is 
this second trajectory merely a continuation of the 
first or is a different particle being observed? • 

Using t to represent time for the “first” particle 
(so that at 1 q, / = 0), its equations of motion are 

x(t) = Ky cos Xyt + K 2 sin Aj t 

+ K 3 cos X 2 t + K 4 sin X 2 t (7a) 

y(t) = My (K 2 cos A^r — Ky sin Xjf^ 

+ M 2 (k 4 cos X 2 r — K 2 sin X 2 r) (7b) 


Using r to represent time for the “second” particle 
(so that at f l5 x = 0), its equations of motion are 

x(x) = ATj cos A^x + K * 2 sin AjX 

+ K 2 cos X 2 t + K ^ sin X 2 r (8a) 

y(r) = My(K 2 cos AjX — K* sin Xjx) 

+ M 2 (a^ cos X 2 x — sin X 2 x^ (8b) 

The relation between rand x is t = x + ty. 

Equations (7a) and (7b) for the original trajectory 
[x(r),y(r)] when viewed in the time frame x become 

x'(t + = Afj cosXj^r + + A' 2 sinXj^T + 

+ KyCosXy^r + t ^ + A' 4 sinX 2 ^r + 

(9a) 

y(r + cosXj^x + - A'j sinXj^x + tjjiJ 

+ My [A 4 cosX 2 ^x + - A ' 3 sin X 2 ^x + 

(9b) 

Using the addition formulas for sin X(x -f tj) and cos 
\(t + fj) in equations (9 a) and (9b) and equating the 
coefficients of cos Xjx, sin Xjx, cos X 2 x, and sin X 2 x 
to the corresponding quantities in equations (8a) and 
(8b), 2 it is found that 

A!^ = A^cosA^ + K 2 sin Xyty (10a) 

A* = K 2 cosXyty — Ky sin A^Xj (10b) 

A| = cos X 2 rj + K 4 sin X 2 ty (10c) 

A* = K a cos X 2 ty - K 3 sin X 2 ty (lOd) 

2 A sufficient condition to permit the equating of the corre- 
sponding coefficients is that X t and X 2 be incommensurate. This is 
assumed here and throughout the remainder of this paper. 
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However, equation (10) implies that 


following relations: 


A* 2 + A* 2 = Aj 2 + A 2 2 (11a) 

, *1 

COS 0J - — — 
J 1 

sin 0j 

*2 
" h 

and 

K 1 

COS 0 4 - J. 

sin 0* 

tr|v 

11 

K* 2 + A* 2 = A 3 2 + A 4 2 (lib) 

. *3 

COS <p 0 = -7— 

J 2 

sin 0 2 

= ^4_ 
J 2 

Therefore, Aj 2 + K 2 2 and A 3 2 + K 4 2 are invariants of 
the trajectory. The notations J^ 2 = Aj 2 + A 2 2 and J 2 2 
= A 3 2 + A 4 2 shall be used. 

* *3 

cos h = 7T 
J 2 

sin 0* 

_ a 4 

j* 

J 2 


4.2 The Third Invariant 

It has been shown that the quantities 7j and J 2 are 
invariants of the trajectory. The question of addi- 
tional invariants is pursued. 

From equations (10a) to (lOd), it is found that 


(Remember that J* = J\ and J 2 = J 2 . Also, only 0* 
and 0 2 are functions of fj.) Using this notation. 


cos^r^ = cos 0^ cos 0* + sin 0^ sin 0* 


*1*1 + K 2 K\ 

cos Vl = 71 

J \ 


If If T' ^ 

K 2 K l K \ K 2 

sin A.,/, = 

11 J 2 

■'l 


= cos(0j - 0*) 

sin^X^j^ = sin0jcos0* — cos 03 sin 0 i 
= sin ^0 2 - 0*) 


cos X 2 r 1 


sin X 2 fj 


k 3 k* + a 4 a* 

J 2 
J 2 

A4A3 A 3 A 4 

J 2 
J 2 


Therefore, there exists an integer n such that 

X l^l = _ 4>i + 2mr (12a) 

It is claimed that n is zero. Consider 


2im = X^ — 0j + 0* (12b) 

(Note that, on the right sides of these equations, only 
AT*, K 2 , A 3 * and A 4 * depend on t^\ the other symbols 
do not.) 

Four angles 0 b 0*, 0 2 , and 0 2 are defined by the The symbol 0j, which depends on the “initial” con- 
ditions at /q, is independent of fj and for our purpose 
is constant. The symbol 0* depends on the “initial” 
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conditions at <f>* is a continuous function of t x . 
Furthermore, if t x = 0, = <f> x . Therefore, the right 

side of equation (12b) is zero for t x = 0. Thus, 
n = 0 when t x = 0. But n is an integer and the right 
side of equation (12b) is a continuous function of t x . 
This implies that n = 0 and equation (12a) becomes 


Vi = h ~ (12c) 

A similar argument shows that 


Vl = <t> 2 ~ $*2 ( 12d ) 

Eliminating t x , 

Vl - X^ 2 = - X^* (13) 


= m x k 2 + m 2 k 4 

(14c) 

= -XjMjAj - X 2 A/ 2 A 3 

(14d) 


When these values are substituted into equation (3c), 
that equation becomes 

(Xj 2 - “2 2 V)*2 2 + 0l 2 V - ^1 2 K 2 

+ ( X 2 2 - “2 2 ^2 2 K 2 + ( W - “l 2 K 2 

+ 2 ^XjX 2 - co 2 2 M l M 7 )K 2 K 4 
+ 2 (X 1 X 2 M 1 M 2 - w x 2 ^K x K 2 
= -C (15) 


The symbol used shall be J 3 = X 2 <f> x — X x (f) 2 . 

Equations (11) and (13) show that if the same tra- 
jectory is observed at two different times, the quan- 
tities 7j, J 2 , and / 3 are invariants. Conversely, if two 
trajectories are given for which the J x , J 2 , and / 3 of 
the first are equal, respectively, to the /*, J 2 , and / 3 of 
the second, then the two trajectories are not distinct; 
one is a continuation at some later time of the other. 
Furthermore, this later time can be determined from 
either equation (12c) or equation (12d). 


However, since M X M 2 = co 1 /ta 2 and XjX 2 = tojcu 2 , 
the coefficients of K 2 K 4 and K X K 2 are zero. 

The coefficients of K x 2 and K 2 2 are shown to be 
equal by the following. From section 3.4, 


co 1 2 + Xj 2 2X x 


4.3 On the Relation Between J 2 , and C, the so that 
Jacobi Constant 


The Jacobi integral, equation (3c), is valid for all 
values of the time t and in particular for t = 0. The 
following relationships between x 0 , y 0 , x 0 , and y 0 and 
K x , K 2 , A 3 , and K 4 can be found from equations (4a) 
and (4b) and their derivatives with respect to time. 



co 1 2 + Xj 2 
co 2 2 + Xj 2 


<x> 2 M^ + XjVj 2 = co 1 2 + X x 2 


X o = K l + K 3 


X o ^1^2 + ^2^4 


(14a) 

(14b) 


\ 2m x 2 ~ <^i 2 = Xj 2 - u 2 M 2 

A similar calculation shows that the coefficients of 
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K 2 2 and K 4 2 are equal. Thus, equation (15) can be 
written as 

(W - X 1 2 K 2 + (“ 2 V - h % 2 = C 

(16) 

5.0 ILLUSTRATED TRAJECTORIES 

Trajectories associated with different initial condi- 
tions are explored, in a systematic way, in this sec- 
tion. At a first level of classification, trajectories were 
differentiated on the basis of the value of the Jacobi 
constant C. In sections 5.1, 5.2, and 5.3, trajectories 
for which C = — 1 , C = +1 , and C = 0, respectively, 
are examined. 


5.1 Trajectories for Which C — — 1 

Included among the trajectories for which C = 
— 1 are the rest trajectories, those for which the start- 
ing position is on the Jacobi ellipse. 

t , 2 „ 2 , 2 , 2 _ , 

W 1 x o + y o ~ 1 


Typical trajectories are drawn for y = 0°, 15°, 
30°, . . . , 345° in figure 1. The value of n = 1/82.45 
was chosen, from which Xj, X 2 , w 2 , A/j, and Af 2 

were calculated. 

For the specific initial conditions, Aj, K 2 , K 2 , K 4 , 
and /j, J 2 were also calculated. Each graph is scaled 
so that each mark represents 5 units. The innermost 
curve (in black) is the Jacobi ellipse; its semi-major 
axis equals 1/coj and its semi-minor axis equals l/to 2 . 
The trajectory is drawn in red for 0 ^ t ^ 50. 3 The 
blue curve will be discussed in section 6. 


5.2 Trajectories for Which C — +1 

Among the trajectories for which C = +1 are 
those which start, with velocity different from zero, 
at the libration point x = 0, y = 0. Under these con- 
ditions, the Jacobi integral becomes 

x 2 + y 2 = l 

o 


These trajectories constitute a one-parameter family 
of solutions that can be parameterized by x 0 = cos y 
and y 0 — sin y, where y is the parameter. The values 
of Kj are 


These constitute a one-parameter family of solutions 
that can be parameterized by x 0 = (1/wj) cos y and 
y 0 = (l/o> 2 ) sin y, where y is the parameter. The 
values of K, are 


K, = 


(to. 2 + X 9 2 ) cos y 

'' «i(s 2 -V) 

2X 2 sin y 


K 2 = 


>l( ? 


co,(X 2 z - X x 2 ) 


*3 = 


*4 =- 


+ Xj 2 ^ cos 7 

' ",(^ 2 2 - V) 

2Xj sin 7 

w l( X 2 2 - X l 2 ) 


= 2 sin 7 

‘ " (* 2 2 - V) 

Xi(cOi 2 + X 2 2 )cos7 
^ = " -1 2 ( X 2 2 - X l 2 ) 

K = 2 s^ 7 

3__ ( X 2 2 - X t 2 ) 

_ X 2( W 1 2 + X l 2 ) cos T 
4_ -1 2 ( X 2 2 - X 1 2 ) 


Typical trajectories are drawn for y = 0°, 45°, 90°, 
135°, 180°, 225°, 270°, and 315° in figure 2. 


3 This corresponds to approximately 217.422 days. 
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5.3 Trajectories for Which C — 0 

The third possible value for Cis zero. Under these 
conditions, the Jacobi integral becomes 


y 2 l 2 — . . 2 2 t 2-, 2 


The trajectories to be used as models will be those 
that start on the ellipse 


“ i 2 ^ 2 + <* 2 2 y 0 2 = “ 2 


with the proper initial velocity dictated by the equa- 
tion 


V + V = a 


These conditions may be parameterized by x 0 = 
(a/ojj) cos /3, y 0 = (a/w 2 ) sin /J, x 0 = a cos y , and j> 0 
= a sin y. The values of K, are 


a = 1, most of the trajectories lay within the bounds 
of the size of the paper available (figs. 3(a) to 3(h)). 
For those trajectories that exceeded these bounds, 
the drawings were rescaled by choosing a = 0.5; 
these trajectories are given in green (figs. 3(i) to 
3(n)). 


6.0 ENVELOPES 

For the trajectories drawn, there is an exterior 
bounding curve for each. In section 6.1, a geometric 
description is given that enables determination of 
this bounding curve. The equation for the bounding 
curve, found as an envelope, is given in section 6.2. 
In section 6.3, orbits of spaceship observatories are 
described. An approximation to the bounding curve 
had been known; in section 6.4, this approximation 
is compared to the envelope determined in section 
6.2. From the analysis given in section 6.2, an inner 
region from which the trajectory may be excluded 
can be determined; this is done in section 6.5. It was 
observed that, for all the trajectories drawn, J 2 ^ J\. 
The question thus arises: Is this always true or are 
there initial conditions for which J\ > J 2 1 This is 
answered affirmatively in section 6.6 and trajectories 
for a select set of initial conditions are drawn. 


6.1 A Geometric Description 


aj^Wj 2 + X 2 2 )cos/ 3 + 2cOj sinvj 

w l( X 2 2 - X l 2 ) 

aXj 2X 2 2 sin0 - (a>j 2 + X 2 2 )w 2 cos7 

2 2 2 

co i co 2 ^2 ^1 

aJjtOj 2 + Xj 2 ) cos/3 + 2COJ sin 7 J 

w l( X 2 2 - X l 2 ) 

aX 2 (^l 2 + X i w 2 cos ^ — 2X 1 2 sin|3 

^2 X 2 


Typical trajectories are shown in figure 3 for /3,-y = 
0°, 45°, 90°, 135°, 180°, 225°, 270°, and 315°. With 


Consider a ray emanating from the libration point 
and the intersections of the trajectory with it. If the 
ray is referenced by an angle 0, let r(0) be the least 
upper bound of all the intersections. Then, r = rifi) is 
the bounding curve for the trajectory. 

On a fixed ray r, the distances from the libration 
point x = y = 0 to successive intersections with the 
trajectory seem randomly generated. To proceed 
directly from this sequence of distances to their least 
upper bound seems hopeless; an alternate descrip- 
tion is needed. 

Consider a particle P moving according to equa- 
tions (4a) and (4b) where the K,- are given by equa- 
tions (5a) to (5d) for the initial conditions (x 0 ,y 0 ), 
(x 0 , y 0 ). At t = 0, when Fis launched, a spaceship Sj 
is also launched. In describing the motion of Sj, 
whose equations of motion will also be equations 
(4a) and (4b), K*, K 2 , K 2 , and shall be used to 
denote the coefficients. The initial conditions for S) 
will be represented by (x 0 * y*), (Xo,y* 0 ). The space- 


10 



ship Sj is launched so that 


As viewed by an observer on S b the particle P de- 
scribes the ellipse £3 


y 


* 

O 



x = A 3 cos \ 2 t + K 4 sin X 2 r 
= J 2 cos(X 2 r — (p 2 ) 


and 



i.e., the initial conditions for spaceship are chosen 
so as to make £ 3 *= A 4 * = 0. 

Under these conditions, the trajectory of the 
spaceship Sj is given by 


y = M 2 (K 4 cosXjt — £3 sin X 2 r) 
= M 2 / 2 sin(X 2 r - 0 2 ) 


The semi-major and semi-minor axes for £j are Jy 
and A/j/j, respectively; for £ 3 , they are J 2 and M 2 J 2 , 
respectively. An observer fixed in the (x,y) system 
and looking at .Sj and P can visualize the ellipse £3 
being attached to Sj and moving with it. The center 
of £3 is at £j and the axes are parallel to the axes of 
£ b as shown in the following sketch. 


x* = K* cos Xj t + K 2 sin X^t 


The ellipse E 3 rides on the ellipse £j 


y* = M ^ (K^ cos X^t - sin X^ 


where A* = x* and K 2 = x*/X [ . If x* and x* are 
chosen as x* 0 = Aj and x 0 = Xj A 2 , then Aj = Aj 
and A 2 = A 2 . Thus, for appropriately chosen initial 
conditions, the trajectory of spaceship Sj is the 
ellipse £j 


x* = A x cos Xj? + K 2 sin Xjt 


*t 



x 


= J 1 cos (Xjf - 0j) 
y* = (A 2 cos Xjf - Aj sin Xjt) 

= sin^Xjf — <p ^ 


Again, consider the particle £ moving according to 
equations (4a) and (4b) where the A,- are given by 
equations (5a) to (5d) for initial conditions (x 0 ,y 0 ), 
(*(kVo)- At t = 0, when Pis launched, a second space- 
ship S 2 is also launched. In describing the motion of 
S 2 , whose equations of motion shall also be equations 
(4a) and (4b), A[*, K* 2 \ and K* 4 * shall be used 
to denote the coefficients. The initial conditions for 
S 2 will be represented by (x^*, yS*), (xj*, y* 0 *). The 
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spaceship S 2 is launched so that 



i.e., the initial conditions for spaceship Sj are chosen 
so as to make K^* = K^* = 0. The trajectory for 
spaceship 5^ is given by the equations 


x** = A"**cos\ 2 r + A**sinX 2 r 
>-** = M 2 (Kl*cos\ 2 t - A**sinX 2 r) 


where K* 3 * = x* 0 * and K* A * = x* 0 */\ 2 . If x* 0 * and x* 0 * 
are chosen as x%* = K 3 and = X 2 Af 4 , then K$* = 
K 3 and K^* = K 4 . Thus, for appropriately chosen 
initial conditions, the trajectory of spaceship S 2 is the 
ellipse E 2 

x ** = K 3 cos X 2 r + K 4 sin X 2 r 
= J 2 cos (\ 2 t — 4> 2 ^ 
y** = M 2 (K 4 cos X 2 r — sin X 2 t) 

= M 2 J 2 sin (X 2 t - 0 2 ) 


and, when viewed by an observer on S 2 , the particle 


P describes the ellipse £ 4 


x = cos Xjt + K 2 sin Xjt 
= J 3 cos (X^t - 0j) 

y = M 3 (K 2 cos Xjt — sinXjt) 
= MyJ^ sin(Xjt — 


The semi-major and semi-minor axes for E 2 are J 2 
and M 2 J 2 , respectively; for £ 4 , they are and M X J X , 
respectively. The geometry for this situation is il- 
lustrated in the following sketch. 


The ellipse E 4 rides on the ellipse E 2 


y 



In the subsequent discussions, the description will 
be based on the image of the ellipse with the smaller 
semi-major axis riding on that spaceship whose tra- 
jectory has the larger semi-major axis. In effect, J x 
and J 2 are being compared, since M x = 0.49 . . . and 
M 2 = 0.19 . . . for fi = 1/82.45. 
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6.2 The Equation of the Envelope 

When a fixed ray emanating from the origin inter- 
sects the smaller ellipse as it is being carried pig- 
gyback on the larger ellipse, the farthest point of in- 
tersection is an upper bound for all intersections of 
that ray and the trajectory. This farthest point of in- 
tersection can be found quantitatively. 

Assume that > J 2 . The equation of the ellipse 
which the spaceship £j traverses is x = J x cos a and 
y = M X J X sin a, where a, the eccentric angle, is con- 
sidered a parameter. 4 The equation of the smaller 
ellipse on whose perimeter the particle rides is 


(x - /j cosa^ 2 - Mj/j sina^ 


M 2 / 2 
m 2 J 2 


= 1 


(17) 


If the envelope of this one-parameter family of 
ellipses is determined, a bound for the trajectory will 
be found. To find the envelope, first differentiate 
equation (17) with respect to a. This produces 


When these two values of y are substituted into equa- 
tion (18), two values of x are determined. 

M.J 0 cos a 

x # = J 1 cos a ± ■ - (19b) 

-Wa/j 2 cos 2 a + sin 2 a 


It is important to remember that the plus signs go 
together as do the minus signs. Equations (19a) and 
(19b) are the parametric equations for the envelope. 
Using the plus signs, the bounding curve is drawn in 
blue in figures 1 to 3. 

The analysis for J 2 5 J x parallels that given above 
and leads to the equations 5 


M~J, cos a 
2 1 

= J 2 cos a ± ~ . — 

'\M 2 2 cos 2 a + M^ 2 sin 2 a 


y # = M^J 2 sin a ± 


Mj 2 /j sin a 


■^M 2 2 cos 2 a + Mj 2 sin 2 a 


— -(x - J j cos a) sin a 




(y 


A/j/j sin a) cos a = 0 


6.3 The Orbit of the Spaceship 


which can be simplified to 


(* 


J j cos a) sin a 


M , 


A/, 


-(y - sin a) 


cos a 


(18) 


Each spaceship, because of the special way the ini- 
tial conditions were chosen (i.e., for Sj, J 2 = 0; for 
Sj, 7j = 0), has a periodic orbit, an ellipse. These 
elliptic orbits are described in reference 2 (pp. 
258-261). 


6.4 An Approximation to the Envelope 


It remains to solve equations (17) and (18) for x and 
y. Solve equation (18) for (x — J\ cos a) and 
substitute this into equation (17). The resulting equa- 
tion can be solved for y . ; there are two solutions: 

M^ 2 J 0 sin a 

y # = M l J 1 sin a ± . (19a) 

yjM i cos l a + A/ 2 sin z a 


An approximation to the envelope developed in 
section 6.2 can be obtained by arbitrarily assuming 
the envelope to be an ellipse with semi-major axis 
+ J 2 and semi-minor axis M ^ + M 2 J 2 (ref. 2, p. 
283). A comparison of this approximate bound and 
the bound found in section 6.2 is shown in figure 4. 
The bounding curve is drawn in blue and the ellipse 
in red. At the scale in which the figures are drawn, 
the approximation is remarkable. 


4 Note that a = Xj/ — <f> y 


5 Note that, in this case, a = \ 2 t — <t> 2 . 
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6.5 An Inner Bounding Curve 

In section 6.2, an outer bounding curve was ob- 
tained when the plus sign was used in equations 
(19a) and (19b). From the geometry of a family of 
ellipses whose centers lie on a Fixed ellipse, one 
suspects that, under the right conditions, there may 
be an inner region from which the trajectory is ex- 
cluded. The existence of such an inner region de- 
pends on the relative magnitude of J x and J 2 and the 
relative magnitude of M X J X and M 2 J 2 . 

The case of C = 0 is disposed of easily. From 
equation (16), it is found that 



and for pt = 1/82.45, 


J x = 0.3515131989732791 A, 

and 


(Myf^ ) = 0.8 87827 1 3962 1 7 024 

In this case (refer to the sketch on page 12), there 
is always an inner region from which the trajectory is 
excluded. Equations (19a) and (19b), when the 
minus sign is used, produce a curve that looks like 
the one shown in the following sketch. 


gular regions, AOjD and B0 3 C, are notpart of the ex- 
cluded region; the excluded region is bounded by 
O 1 O 2 O 3 O 4 . A typical situation is shown in Figure 5. 
In Figure 6 , the family of smaller ellipses is drawn in 
black, the trajectory in green, and the bounding 
curves in red. 

For C = +1, there may be no inner excluded 
region, as can be seen in Figure 7. The inner envelope 
in this illustration, drawn in green, does not encom- 
pass a region from which the trajectory is excluded. 
In the next section, initial conditions are given for 
which there does exist an inner region from which 
the trajectory is excluded. 

For C = —1, there will always be an inner ex- 
cluded region that includes the Jacobi ellipse. 


6.6 On the Relative Magnitude of /, and J 2 

Equation (16) can be rewritten as 



For fx = 1/82.45, 


7j 2 = (0.1235615290524281)7 2 2 + (5.1 82 1 89 1 801 03332)C 


For all the Figures drawn thus far, J x ^ J 2 ; are there 
initial conditions for which J x > J{1 If C = 0, 


J x = (0.3515131989732791)^2 



For the parameter a (a going from 0 to 27r), 
the sequence to be followed in tracing the curve is 
QA 0 1 020 3 BC 0 3 0 4 0 1 DQ. However, the two trian- 


and clearly J x < J 2 . Similarly, if C = — 1, 


J x < (0.3515131989732791)^ 2 

Thus, if any trajectories are to be found for which 
J x > / 2 , they can only be among those for which 
C = +1. 

It is not intended to exhibit samples of all trajecto- 
ries for which J x > J 2 , only some of a special set. If J x 
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> J 2 , then > J 2 and + AT 2 2 > A^ 2 + AT 4 2 . 
Using the parametric representation of equation (22) 
in the appendix, 


(coj 2 + X 2 2 )cosha sinh |3 + 2io^ cosh a cosh (3 sin 7 
“l( X 2 2 - X l 2 ) 

Xj |"2X 2 2 sinha — co^cjj 2 + X 2 2 ) cosh a cosh (3 cos 


2 2 
w, 


C ° 2 ( X 2 2 - X l 2 ) 


^co j 2 + Xj 2 ^cosh a sinh (3 + 2cj^ cosh a cosh (3 sin 7 
co ^^ 2 - Xj 2 ) 

X 2 ["wj^coj 2 + Xj^cosh a cosh (3 cos 7 - 2Xj 2 sinha 


4 2 

U| 


W 2 ( X 2 2 - X l 2 ) 


Special values are chosen for a and y; a = 0,y 
90°. For these values, 


(o}y 2 + X 2 2 )sinh |3 + 2 oj^ cosh /3 


"l(*2 2 - X l 2 ) 


k 2 = 0 


+ Aj 2 ^sinh (3 + 2 cjj cosh (3 

A ' 3 ‘ Mh 2 - \ 2 ) 


K t = 0 


Thus, > J 2 implies that | | > | K 3 \ . 


Case 1 


[(V + A 2 2 ^sinh (3 + 2 wj cosh/ij > 

— j 2 + Xj 2 ^sinh /3 + 2 cjj cosh/ij 


This implies that 


tanh j 3 > — 


2cOj 2 + Xj 2 + A 2 2 


or that /3 > -0.7343384813667741. 


Case 2 


I sinh /3 + 2 co, 


This implies that 


(A 2 2 - Aj 2 )sinh /3 > 0 


or, since (A 2 2 — A t 2 ) < 0, that sinh /3 < 0. However, 
sinh (3 < 0 for /3 < 0. Therefore, we find that J x > J 2 
for a = 0, -0.73 < /3 < 0, y = 90°. Trajectories for 
/3 = —0.1 through —0.7 are given in figure 8. The 
shape of the interior excluded region in figure 8, as 
shown in the following sketch, is similar to that of 
the sketch shown in section 6.5 rotated 90°. 
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7.0 CONCLUDING REMARKS 

The derivation of the finite equations of motion 
from the differential equations that describe the mo- 
tion of a point particle in the neighborhood of a 
triangular libration point ( L 4 or L s ) in the linearized, 
restricted problem of three bodies is well known. 
From these equations, three integral invariants of the 
motion were derived; the Jacobi integral is expressed 
linearly in terms of two of these. 

The trajectories for varied initial conditions were 
drawn. The motion of a point particle as it traversed 
its trajectory was described in elementary geometric 
terms. 

Each trajectory has an exterior boundary curve; 
its equation was found. An approximation to this 
boundary curve had been known; it is a well-defined 
ellipse. For different initial conditions, the boundary 


curve and the ellipse were compared graphically; the 
agreement seems remarkable. For certain initial con- 
ditions, there is an interior region from which the 
trajectory is excluded; the equation of the boundary 
of this region was found. 

Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, September 28, 1979 
953-36-00-00-72 
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Appendix 

Parameterizations for the Jacobi Ellipse 


A convenient method for choosing the initial 
values of the position and velocity of the particle can 
be obtained from the Jacobi integral: 

V + ^ 0 2 - W - ^o 2 = c (20) 


If C — — 1 , equation (20) can be parameterized by 


cosh a cosh B cos y 
° oj 1 ' 

(21a) 

v = — — cosh a cosh B sin y 
° w 2 

(21b) 

x Q = cosh a sinh /3 

(21c) 

y Q = sinh a 

(2 Id) 


For the rest trajectories (sec. 5.1), a = /3 = 0, y ar- 
bitrary. If the initial velocity is not zero, a and /3 can 
be determined from equations (21c) and (21d) and y 
from equations (21a) and (21b). 

If C = +1, equation (20) can be parameterized by 


For trajectories that start at the libration point (x 0 = 
y 0 ~ 0), a = (3 = 0, y arbitrary. If the initial position 
is different from the origin, a and /3 can be deter- 
mined from equations (22a) and (22b) and y from 
equations (22c) and (22d). 

If C = 0, equation (20) can be parameterized by 


x n = — a cos B 

° CJj 

(23a) 

y a = 7^-“sin/3 

(23b) 

x Q = a cos y 

(23c) 

y o = a cos y 

(23d) 

Included among these trajectories is the solution cor- 
responding to a = 0; i.e., the particle is at the libra- 
tion point with zero initial velocity and remains there 

indefinitely. Leaving this solution 
parameterization is as. given above. 

aside, the 


= — — cosh a sinh B 
W 1 

(22a) 

= — sinh a 
w 2 

(22b) 

= cosh a cosh j5 cos y 

(22c) 

= cosh a cosh /? sin y 

(22d) 
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+ 




x 0 = 6.062834463 
x o = 0 

Ki = -0.8547469351 
K 3 = 6.91 7581398 
J. = 0.8547469351 


y 0 -o 
y 0 = o 
K 2 = 0 
K< = 0 

Js = 6.917581398 


(a) a = 0, 0 = 0, y = 0°. 


FIGURE 1. — Typical trajectories for which C= — 1. 



x 0 = 5.856248389 
x o = 0 

K, = -0.8256221396 
Ka = 6.681 870528 
Jt = 1.404920191 


y 0 = 

y 0 - 

k 2 
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